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The treatise seeks to find the optimum shape of a thin 
Permeunder the Influence of lateral buckling. The specific 
formulation is by means of the isoperimetric problem of the 
calculus of variations. An energy method approach yields 
the governing equations for the problem. An analytic 
solution was not obtained due to the nonlinearity and 
coupling conditions of the equations. 

A variable height lateral buckling problem is formulated 
as an alternate attempt to find the optimum design. Through 
Bessel equation analysis numerous designs are obtained and 
the resulting buckling loads are calculated. The largest 
buckling load corresponds to a beam design which appears to 


be very close to the cptimum shape. 





TABLE OF CONTENTS 


ir. ete CO fe ee SE eee 9 
A. THE OPTIMIZATION PROBLEM ------ GS le 9 

Eee Cor eOr i Tiaels —--—----- =SSacsee == 13 

i sor Latie OF OPTIMAL STRUCTURAL DESDGN ------------ 13 
A. THE BASIC PROBLEM OF STRUCTURAL MECHANICS --- 13 

ee) ee Od ee 14 
CoO a nih Oe ReOEmE = —--—=—— 15 

D. POTENTIAL ENERGY ---------------------------- 16 

fe Pine irrEOneMimmenun POTMINTT AL —------—----—— a 

Po oinetnAl Ortiii7ATkON PROBLEM ——————————— 18 
nh Aten, Bue [NGwPRO sli —---s===——------ = 22 
a CHER AL) DESChhPnieah ——————--———-————— ——_-= ___==— 22 

B. DEVELOPMENT OF CLASSICAL SOLUTION ----------- 23 

IV. A VARIABLE HEIGHT BUCKLING PROBLEM -------------- 32 
Pe DESCRIPTION OP 2Hh PROREBM ------------—-__-- 32 

1. Non-integer Value v --------------------- 35 

2. Integer Value of v ---------------------- So) 

B. CALCULATION OF AMMENDED LOADS --------------- 38 

C. EXTENSION OF THE PROBLEM -------------------- 43 

D. LATERAL BUCKLING OF A CANTELEVER ------------ 45 

E. COMPARATIVE ANALYSIS OF SOLUTION ------------ 47 

F. VALUES OF n GREATER THAN ONE ---------------- 52 

ye Selmemieneis —-----__--__-_____-_- 54 
APPENDIX A: A STRUCTURAL OPTIMIZATION PROBLEM -------- 56 





Poe Bs EE PPECT OF THE LOAD P AT A DISTANCE 


Hee er E@sTeROmD —-----=-----L=.---..-____ 62 
APPENDIX C: TRANSFORMATION OF EQUATION (4.11) 

 aene miEenmuIeN ---—---.—-.....___....- 6 
Pere ee: ShcSEL EQUATION CURVES --<---------------- 66 
ieee Eb: MVARDABLE HEELGHT DESIGNS --------~-~-----~-- 76 
Corwen eRe eR mms =-=--5----------------=..----------..- 108 
(ii eo Ben eGR S) -o--2->—---..-....-..._--...._.---~- 123 
ee pier rEyEmenwenmeT —--------------_-.~--~----_-~_ 125 
et en ee | ee ~~ 126 





1. Column Shapes Analyzed ---------------------------- 18 
©. Buckled Beam ------------------- Hn en ee 22 
3. Beam Loaded at the Centroid ----------------------- 23 
Lateral Buckling Configuration +------------------- 25 
5. Variable Height Beam -----------------—---------+----- 29 
6. Variable Height Beam -------------~-------—---------- 32 
~. Convergence of Power Law -~-------------<----------- 39 
8. Variable Height Beam (n=0.14) --------------------- oO 
9. Graph P* vs n (S.S. beam - case 1) ---------------- 43 
10. Variable Height Beam ------------------------------ 4 yy 
11. Graph P* vs n (S.S. beam — case 2) -----~----------- 46 
12. Graph P* vs n (cantilever beam ) ------------------ 48 
13. Variable Height Cantilever ------------------------ LQ 
14. Moment Configuration ------------------------------ 51 
15. Variable Width Beam ------------------------------- 56 
16. Optimum Shape Variable Width Design --------------- 61 
17. Effect of P at the Top of the Beam ----------~----- 63 
18.-27. Bessel Equation Curves ------------------------ 66-75 
28.-59. Variable Height Beam Designs ------------------ 76-107 


Lis# OF FGURES 


ope a fF 


EtSiSOr SYMBOLS 


Cress secturemal area 
€Cerewant of integration 
variable width 

constant width 

maveriral constant 
maeerial constant 
domain 

Youngs "s’modulus 
fumetional 

filme t von 

shear modulus - 
variable height 
constant height 
paeearaletunction 
WESFOr SubsSeript 

polar moment of inertia 
Vector Subscript 
dimensionless constant 
dependent variable subscript 
length 

moment 

column load 

exponent for design 
COnGemutraved load 


Unmeeroarmn load 





Ss surface 


ah total potential energy 

U strain energy 

V potential energy external forces 
LA volume 


Leng displacements 


ey; 2 eCoordmmave axes 

a constant of integration 

B Bessel equation parameter 
Y (ij ORL) 

6 variation 

€ strain 


Boat) g principle axes 


r . bagreanewanemuiltiplier - 
v Bessel function order 

fe) Bessel equation argument 
O Swieeos 

y Sect lie Strain energy 





ACKNOWLEDGEMENTS 


The author wishes to express his appreciation to his 
Geer eosoctave Professor David Salinas for his construc- 
tive supervision and enthusiastic encouragement throughout 
this work. He also wishes to thank Associate Professor 
Craig Comstock for his help in overcoming mathematical 
difficulties encountered during the thesis. 

The Naval Postgraduate School Computer Center provided 


maciiities for computer work. 





1. DNERODUCTION 


Pee OPTEMEZATION PROBLEM 

The performance of an elastic body is measured by the 
maximum load which the body is able to withstand. The 
weer Zav1on problem considered in this treatise is that 
Pane orcrmnine the specific design of an elastic structure 
such that the structure will exhibit a maximum performance 
capacity. Herein, the term design refers to the actual 
paystca l ase purton of material within an elastic body. 

mae Optimization problem is formulated by the calculus 
of variations. This method offers certain advantages in 
the development of a theory for optimization of an elastic 
boay. Phe specific formulation of the problem is by means 
of the isoperimetric problem of the calculus of variations. 
By this method the problem becomes an extremization of The 
total potential energy of the system with a superimposed 
Vomeume Constraint. The extremization of the funetional, 
augmented by Lagrangian multipliers, leads to a set of 
Euler equations which govern the design of the structure. 

The calculus of variations is no new technique for 
obtaining an optimum design. Lagrange [1] formulated an 
opramum column problem by this means in 1773. He sought to 
determine the shape of a column for which the buckling load 
was largest among all columns of given length and volume. 


Although he failed in obtaining the correct solution, his 





recognition of the problem belonging to the calculus of 
variations inspired others to continue research in this 
field. Clausen [2] later provided a soiution to this prob- 
fem or columms of circular cross sections. Keller [3] 
mecent liye shewed that the strongest column #s not ome ofa 
circular cross section, but instead, is one which has an 
Gul mewereal triangle as a cross section. 

The isoperimetric problem for optimum design received 
an additional stimulus from Keller and Tadjbakbsh [4]. This 
work provided a variational technique aoe 7 oneal Volume, 
elastic design of columns. Taylor [5] treated the same 
problem by means of an energy method approach. In his work 
Taylor showed that the same governing equations may be 
Cvereumcd by direct extremization of a funetional, which 
consists of the system potential energy and a volume con- 
Sweeny.  Layler™Ss approach is the esseftial- method used 
for optimization in this thesis. The method is particularly 
useful since it provides a simple and direct means for iso- 
perimetric problem formulation. 

An extension of the Taylor approach was accomplished a 
Salimes [OjPand provides some conditions for the validity 
of the potential energy formulation. Additionally, the 
author demonstrates the equivalence of maximum load, minimum 
volume optimum designs and presentsa general formulation for 


Opummezation of elastic structures. 
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Pee cee? OFTHE THESIS 

Wsing the Taylor approach, the optimization problem is 
discussed and developed by use of Lagrangian multipliers 
and the calculus of variations. The resulting set of Euler 
equations are seen to be the system equilibrium equations 
and an optimality condition. Example problems are included 
to illustrate the usefulness of the Taylor method. These 
- include a strength problem for a uniformly loaded beam and 
an optimum design problem for a particular case of column 
puck ling. 

The method is then extended to the case of a thin beam 
subjected to lateral buckling. Under lateral buckling, the 
beam is considered to buckle in a plane perpendicular to 
the vertical plane of loading.- The buckling configuration 
consists of rotation and sidewise bending away from the 
MeriZontal axis of symmetry. The resulting equations 
governing an optimum design are found to be non-linear. [In 
Beem Ol, tme equations are coupled in rotation and bending. 
Pil -aevemovs for an analytic solution proved futile. 

An attempt to find the optimum shape has been made by 
as ume Seemebrical conditions for instability and design 
of the beam. A variable height problem has been formulated 
by means of a method described by Federhofer [7]. Numerous 
designs are obtained by use of Bessel equation analysis and 
@a predicted optimum shape is forseen. 

iirovsnemt the discussion all problems are considered 


only for the elastic range of the material. Additionally, 


ul 





Picwmewerial is considered to be homogenous and isotropic. 
All problems analyzed are restricted to reetangular cross 


sections. 


WZ 





een coemo On OPTIMAL STRUCTURAL DESIGN 


eee es PkolC PROBLEM OF STRUCTURAL MECHANICS 


The classical theory of structural mechanics is concerned 
with the behavior of continuous elastic bodies which are 
subjected to specific loading conditions. In general the 
components which define mechanical behavior are the displace- 
ment vector u, and the stress and strain tensors, O55 and 
ay 

The strain tensor is explicitly defined by the geometry 


of deformation of an elastic body. Hence, the strain- 


displacement relation may be represented by 
e = €(u) : (Zell) 


magma the Constitutive laws of elasticity it is known that 


Seress 1S related to strain, that is 
G = Ge) (222) 


Thus, stress and strain are both representable as functions 
of the displacement. In classical theory any one of the 
three quantities, u, o, and € is sufficient for determining 
the other two. 

The mechanical state of a structure is defined by a 
boundary value problem which is comprised of equilibrium 
eGigiwmeans aid tne associated boundary conditions. Thus, a 
Classucal preplem in structural mechanics consists of using 


the constitutive relations between stress, strain and 


IL 





displacement to derive and solve the governing equilibrium 


“oem itonms and their associated boundary cemditions. 


B. CALCULUS OF VARIATIONS 

Problems concerned with the determination of extreme 
values of integrals whose integrands contain unknown func- 
tions below to that field known as the "calculus of 
Variations" [8]. Such an extremal refers to the maximun, 
Pinte or stationary value of a functional. A functional 
fmomoerea! Valued function of @ Tunction. 


Only the integral functional 


mia) = £ F[f(x))dx (253) 
D 


will be considered herein, where D is the n-dimensional 
domain over which the aod an Sse vector x ranges. The 
function which extremizes I[f(x)] is denoted as f*(x) and 
the associated value of the functional is I* (i.e. I* = 
eax) )) . 


inmeGecmerail preplems in structural mechanics concerning 


optimization will assume the form 


te oeenivga” a" dod! jee. ax 2) 


lint 
af: 
for a one dimensional case. Avone dimensional ease is 


generally characterized as having one dimension, its length, 
much larger than the other dimensions. In the above form x 
represents an independent variable. The quantities u, ¢and 
their respective derivatives represent dependent variables. 


In this thesis u and ¢ represent displacement functions or 


14 





Seave Variables. The nature of the specific problem governs 
the number of dependent and independent variables. 
rerhernime® 2 variation with respect toa functional is 
denoted by the symbol 6. From [8] it can be shown that a 
peaumemeary Value of the functional occurs for 6F = 0. This 
operation yields what is known as a set of Euler equations. 


For example, if equation (2.4) is represented by the par- 


eiecular form 


if oem ju eed” A) dx C25) 
D 


then performing the extremal operation yields three Euler 


equations [9]. 





_ oF 9 {oF 9. {or | _ 
a Zia du «Ox en 7 ee eal : (2.6) 
: XxX 
Eee) I Ce 
6 - 02 Gan - 55 3 | = 6 (2.7) 
2 aF _ 
= O a (238) 


plus a set of natural boundary conditions. 


Co bee. LT SOPERIMETRIC PROBLEM 

ee, Sem@erimetricaproblem ingthe calculus.of variations 
provides for extremizing functionals whose dependent vari- 
ables are constrained by an integral relation. The integral 
constraint is accommodated by constant Lagrangian 
medtipliers. 

As an example, suppose that the maximum (or minimum) 


value of an integral 


BES 





oe F Ge: wea ) dx (2.9) 
R } 


Pmiomee, sub,ect tO the condition that a@hother integral 


ef Gx ust! )dx (2 Aon 
R 


is known to have a constant value. Using the method of 


Lagrange multipliers [8] an augmented functional 


ioe ) Teemu.) + GCxsu,u' fydx (251 
ie 


may be constructed and its free extremum considered. The 
necessary condition for a stationary value is the Euler 


eqviacion 
- = oo = CO (2.12) 
where,F has been replaced by F* = (F+AG). 


De SPOrENTIAL ENERGY 
The total potential energy of a conservative elastic 


system can be expressed as 
Uo 2B 


U represents the total strain energy of the structure, and 

V is the potential energy of external forces associated with 
the displacement of the system under any set of given load- 
ie weonmarerons. The total strain energy of the system is 
the volume integral 


Bie. Seic 


ee. 2) aL! 
5 19 fag0¥ ( ) 


where oe amg ¢€ # are components of the stress and strain 
iJ 


tensors. 


16 





The potential energy of external forces is defined as 
the inner product of the external force vector P and the 


associated surface displacement vector u, that is, 
V=efP-°u4ads 
Ss 
where s is the domain over which P is defined. 


E. PRINCIPLE OF MINIMUM POTENTIAL ENERGY 

A powerful tool in the field of structural mechanics is 
the principle of minimum potential energy. This theorem 
states that the total potential energy of a conservative 
elastic system attains a minimum, with respect to variation 
of kinematically admissible displacement functions, when 
the displacement function satisfies equilibrium. 

Mie total potential energy has previously been defined 
in terms of strain energy and external work. These quanti- 
ties may be expressed as functions of the displacement 
vector u. Thus, the equilibrium equations and associated 
boundary conditions are derivable from a variation with 
memmect tO one Gisplacements of a potential energy functional. 


The variational operation 


6 (U-V) = 0 (27 ee) 
can then be performed for which the resulting Euler equations 
physically represent the system equilibrium equations. The 
optimization method considered in the thesis 1S based on 


Ul wee amportant principle. 


aly 





io 2 Se RUCTURAL OPTIMIZATION PROBLEM 

temeiemcrample of structural optimization, consider the 
Pao weeWieG! —Coraiting the optimum shape for column buckling. 
The problem may be characterized by the statement: "For a 
column of given length and volume of material, determine the 
column shape for which the Euler buckling load is a maximum." 

The column shape is considered to be such that all cross 
sections are similar. The moment of inertia I(x) is 


assumed to be related to area, A(x), that is 
EI(x) = kA™(x) (2 


where E represents Young's modulus and k is a constant 
depending on the section considered. For example, the 
moments of inertia of rectangular sections of variable- 
width-constant height, and variable height-constant width 


ame (Fig. 1), 





Figure 1. Column Shapes Analyzed, (a) variable width - 
constant height, (b) variable height - 
constant width. 
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Eh>b (x) 


EI = 13 =) hei) Ga = fi) eric.) 

mor @ese (a) and 
ex = EBRCx) « 443(x) (n = 3) (2.19) 
yy 12 . me 


fomecase (b). 
mere vne Buller Duckling problem the total potential 
energy is 


Q 


a mt 
710 


= \ereiax (2.20) 


where v represents lateral displacement for the buckled 
Soave. Introducing an isoperimetric condition for a con- 


stant volume of material ie gives 


Q 
tauvoyax- = V C2425) 
O 
0 
The volume constraint is accommodated through the 


Hemmation of the functional T*, defined by 


T? T - aN (Ane) 


where A is a constant Lagrangian multiplier. Combining 
22] eeweeOememcd (2922)-the constrained functional is 
Q 


pee f (K 
0 2 


ee ye = 5 vit = heads (22238) 


The governing equations, in addition to the isoperimetric 
equation, are obtained by performing a variation with 


meomecestoOrveand A to obtain 


re 





tf 
ee = 0: (kAtv") + Nv" = 0 ys) 


§,T* = 0: snkatt ees 


st) 
© 


(2525) 


Equation (2.24) is the equilibrium equation and (2.25) an 
Seametity condition. 

The non-linear equations (2.24), (2.25) and (2.21) define 
the optimum design with respect to constant volume - maximum 
load. In general these are extremely difficult to solve 
in closed form. However, an exact solution can be obtained 


for case (a) where n= 1. These governing equations then 


mecome 
(kAv") = Nv" = 0 (2.26) 
Bikey" sD (anata 
Ne 
' oe dx = Vie (22 2am) 


Fomacion (2.27) implies constant bending curvature, i.e. 





vi" = Y2A/k C2 ee) 


Peet icueame (2.28) imto (2226) and integrating twice gives 


© = N (xo Hoe 5 8 6) (2.29) 
~~ ie Z ali a ° 
Sameer y 
x° 
See van 7 1K Cy + C 3x + Cy) ' (2.30) 


If the column is considered to be simply supported, the 


Pouce, Condtyions for displacement are 


meen = v(t) = 0 (23315) 
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mer which C 


3 - %, Cy = 0, and 


2 
v(x) = Yax7k (A (2.32) 
The moment conditions are 
M(o) = M(2) = 0 (23s) 
from which C, = - 7 C, = 0, and 
atx) = NE BY (2.34) 
-- (5 - >. AS 


The critical load is now determined from the isoperimetric 


GCeomaruLlon 
Q 2 
N /X eo 
for which 
7 ~12kKVo _ See 
N optimum 7 — "  - (2. 36) 


Substituting this expression into (2.34) determines the 


optimum design in terms of known constants, that is 


M 2 
A(x) Se i. Cae 
122 


Note that the Lagrange multiplier does not appear in either 
of the expressions for A(x) or N. This is in distinction 
GO Seremeun optimization problems where A is a representa- 
GiOnN Or Maximum strength and remains in the final analysis 


(see Appendix A). 


au 





iii Tea eae BUCKLING. PROBLEM 


#. GENERAL DESCRIPTION 

Me problem of lateral buckling is important in the 
design of beams which have no lateral support and possess a 
shape such that the beam height is large in comparison to 
the width. Such a configuration may become laterally 
unstable when the loading reaches a critical value. The 
PecWieine failure or “buckle” is a combination of lateral 
bending and twisting away from the horizontal axis (Fig. 2). 
Prior to buckling the load needed to cause this failure may 


be considerably less than the corresponding load that would 





Figure 2. Simply Supported Beam in a Lateral Buckled State. 
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cause bending failure. This is due to a large flexural 

rigidity of the beam in the vertical plane of bending as 
compared to a much smaller lateral bending resistance in 
the plane perpendicular to the applied load. 

Boreucimont Intermagive discussion of the historical 
development of the lateral buckling problem is given by 
Bleich [12] and will not be included here. The methods 
used for development of the problem in this thesis are sim- 


ilar to those used by Sechler [13] and Timoshenko [14]. 


B. DEVELOPMENT OF CLASSICAL SOLUTION 

While the solution to the classical lateral buckling 
problem may be found in a number of texts, its inclusion 
herein is considered necessary for understanding the 
descriptive geometry of the problem and later development 
of the problem of a variable height solution. The example 
to be considered is that of a simply supported beam loaded 


at its centroid by a concentrated load P at the origin 


(Fig. 3). The centroidal loading condition is necessary to 


va 





Yo | 


Pere 3. simply Supported Beam Loaded at the Centroid. 
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Sameiity the problem. Loading at any point other than the 
a erera invroduces an additional work termm which leads to 
BSeme Mathematica difficulties. This condition is described 
in Appendix B. 

It is assumed that the load P remains in the vertical 
myeolane Chroughout bucking. The ends of the beam are so 
Constrained as to prevent rotation in the yz plane. The 
buckled configuration is shown in Figure 4. 

The deflection of the beam is defined at any point x 
from the origin by the components u, v and @. Components 
u and v represent the displacement of the centroid in the 
zZ and y directions while $? is the angle of rotation of the 
Gr@es section in the yzZ plane. Thus, the boundary condi- 


tions for the problem are 


uve) = >(2%72) = 0 (3.1) 


- ut(0) = ¢'(0) = 0 (See 
The fixed coordinate axes are labeled x, y and z. In 


the buckled state the principle axes of symmetry at any 
Secvion mn are defined as €, n, and Z. The usual right 
hand rule convention for vectors is used to describe posi- 
tive moments with respect to the described axes. 

By the usual assumptions of small deflection theory the 


equations for principle moments can be written 


EI. oo» EIv"=M (25a) 


a4 








Figure 4. Lateral Buckling Configuration. 


ae. 





US G ~ woe 

EI, : EI u Mn Gcrecicn. 
de _ A 

GJ qe * GJIo My (3.3c) 


These moments are physically described as follows: 


i) M, is the usual bending moment due to the 
5 load P 


ii) M. is the bending moment associated with 
lateral buckling 


1 ) Me is the twisting moment associated with 
ikatviermal buckling. 


With respect to the fixed x, y and Z axes the above 


moments can be represented [13] as 


Me = Me + Mu! = GJ¢o' (324) 

My = M., 1 M¢ = -ET ut (3.4b) 

eae Moe = ELV" (3.4c) 
Let Cl = GJ and Co = Er: Dri merentitabing, tne first 


Peuavton with respect to x and substituting for u" from the 


second results in a general equation which can be written 





M MM 
eee eC (3.5) 
m2 1 1 bee 


For the problem considered let the deflection at the 
beam center (maximum deflection) in the z direction be 


denoted by u The moments at any distance x are determined 


1° 


as 


Ss 
T 


5p 2 (u,-u) (3 fem 


M = 0 CJrbiey 


Bie 





Mee= P/2 @W/2-x) (src) 


Z 


Pieierenpveti ine ana substituting into equation (3.5) 


pee lds 


eta = 2 Om Cater ) 


Introducing the new variable 8 = (2£/2-x) and letting 
@ 


“ = —- equation (3.7) becomes 
ree 
o" + k°B°d =O (3.8) 


Equation (3.8) is now in the form of a Bessel equation for 


which the general solution has the form 


*% ice. k 
¢6=8 [a,J, (5 RB ) + eae (5 B-)] (3.9) 


2 -% 


where Jy and J are Bessel functions of the first kind of 
u 


al 
-*% 


orders % and -% respectively [14]. For the simply supported 


beam the boundary eonditions with respeet to 9¢ are 


¢= 0 at £~ = 0 Cisey = Le) (3.10a) 

¢'=0O0 at B= k/2 fayenexs= -0) 4 ( rey 
Mme Iirst boundary condition yields ay = 0. Differentiation 
of ¢ with respect to 8B now yields 

do _ 3/2 iey,.2 

= ee egy > . ee 


Applying the second boundary condition shows that buckling 


occurs when 


J ke") = Cri) 
i 3. 


ani 





Prom a2 table of Bessel functions [15] the first root of 


Mell) is 


} 
7 7 1.058508 ° (3. 13) 


Pwipsticvuving for K, the buckling load is found to be 


GE 936 enc. 
le = 
2 


Gis Q 


(35a) 


Equation (3.14) represents the exact solution for the prob- 
lem considered and will be used as a comparative value for 
@ritical loads calculated for variable cross sectional 


designs later. 


C. OPTIMIZATION TECHNIQUE APPLIED TO LATERAL BUCKLING 

ihe energy formulation deseribed in Section Il is now 
msedad to develop the governing equations for an optimum shape 
of a thin beam which undergoes lateral buckling. The 
assumptions made during the preceeding section apply. Addi- 
tional assumptions include neglecting the strain energy due 
cor ‘'wanping"” [14] and elongation of the beam. 

The case considered is that of a variable height-constant 
width beam (Fig. 5). The principle moments of inertia with 


meSspect to the buckled state are 


= eC) <= PS A(x) (aie 
ao 
up 
| ae: mec) 
J = 3 b> h(x) S ACx) 18) 


The strain energy due to lateral bending may be written 


28 








Figure 5. Variable Height, Constant Width Beam. 


ie 
U, = } C,Au"" ax (Saeed) 


where Cy = Ebo/12. The strain energy associated with 


twisting of the structure is 
ae 
CoAb'© ax (3.18) 


* o 
where C., = Gb 0/3. 


As the beam buckles the load P is lowered an amount which 


may be represented as [14] 
Le (Oey co nts > eam (se. 9) 


The potential energy of external forces associated with the 


buckling configuration is determined by 


Q/2 
Ne Pf ONC 2 oe) ola (E20) 
0 


The total potential energy of the system is U - V or 


2 


T= f [C,Au"® + C,Ag'® - Pou" (2/2 - x)}dx. (3.21) 
0 


To formulate a functional by the method described earlier, 


20. 





pa TSoOperimetric constraint is introduced for a given volume 


of material, 


f Adx = ve : Gee) 


The augmented functional may be written 


Lf 2 


a= f [C, Au" 
0 


+ Cag '*~ Pou" (3 — x) —~ XA] ax. (3.23) 


Performing a variation with respect to the variables u, $¢ 


and A yields three Euler equations: 


tt 


are + & 7 
oa. = 0 [2C, Au — Po(s = Coes) 
6 ,T# = 0: 20,,(Ab") + Pull (5 = x) =0 (ARs 
6aHs = 0: Cou"™> + C.6'* - i = 0 (3.26) 
‘ mace 5 


Note that the variation with respect to the state variables 
Peands? Fields equilibrium equations. Variation with 
respect to the design variable A introduces an optimality 


equation. The associated boundary conditions are: 
me Ae = 
Nee OG) = 0 Gye 


ice) = o' (Oo) = 0 (2326) 


Integration of equation (3.24) twice and evaluating at the 


boundary conditions determines 


20, Au" - Po (5 Lane oO Gres) 


Imcerration of (3.25) yields the result 


CAg' + 5 u’ (S Joys = (u,-u) = 0 (2530p 
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where Uy represents the maximum deflection at the beam 
midspan. 

The above equilibrium and optimality equations along 
with the isoperimetric condition (3.22) define the optimum 
shape for the structure. These governing equations are noted 
to be non-linear and an attempt to solve them has been 
Maswecessiul. Unfortunately, these equat@ons resist uncoupl- 
ing, which precludes getting to the isoperimetric condition 
as was accomplished easily in Section III. Additionally, no 
computer techniques were found which could solve systems of 
simultaneous nonlinear equations. 

After considerable efforts to solve the above non-linear 
system proved fruitless, an attempt to formulate a variable 
merent solution rom 4a specified polynomial was undertaken 
miemopes Of predicting or converging to am optimum solution. 
mac metnod used and results obtained are included in the 


following section. 
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IV. A VARIABLE HEIGHT BUCKLING PROBLEM 


ee DEoCRIPTRION.ORSTHE PROBLEM 

The problem of lateral buckling of a wariable height, 
constant width beam is now considered. Again, the problem 
to be investigated is that of a simply supported bean, 
loaded at the midspan centroid. However, the height now 


varies according to the power law 
i a ie) CH) 


where om is the height at the center section of the beam 
wee. 6). “The valve of nm is restricted to lie between the 
mieerval O <n < 1 for reasons which become physically 
Obvious in later analysis of the problem. The idea Mone 


Meine ehas type of formulation was based on the results 


y eC -2X/Q) 





Figure 6. Variable Height Beam Based on the 
Power Law mG 2x72) 


Be 





mecainea@eior a Variable height cantilever problem by 
Federhofer [7]. 

Differentiating equation (3.4a) with respect to x and 
mesing that the previous constants C, and Co are now vari- 


able functions of x, reduces (3.4a) to the form 
t" tat = mt iat 1 
Cio" + Ci" ee iut + Mou" Cie) 
Smescicucing (3.4b) into (4.2) for u" and recalling the 


moment conditions from (3.6) (and their respective deriv- 


Miives) reduces (4.2) to 


Lie 2 
Pa x) 
oN + g- ClO! + Gago = 0 (4.3) 
Al a 2 
Betting 
| y= (1- eX (4.4) 


and using the relations 


_ 1 5 
Cy =e Eb- h(x) 
ag 3 


the lateral and flexural stiffness at the beam center may 


and Ox Thus, the new variable stiff- 


nesses can be represented by 


be designated as C 


n 


= a ao 
220 Ciel AO 
Co = C5 Cie is ) = COY (Gl So) 


Using the chain rule 
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7 aa a CES 
ema noting 
X = > (1-y ) 
2 2 
el OME lou 
> as 4 2 = > a Cie) 
x g dy 
pam larly 
1 1 a6 4 on ad (4.8) 
Cy dx dx Qe vy MY. ; 


Eeoceicution of (4.7) and (4.8) into equation (4.3) yields 


Chin aera” sf deep 
aye Y 4 64c°Cs 


The dimensionless parameter ie may be defined 


etl 
= Ss = Ke ee 16D) } 
64c?c° 


eee 
Equation (4.9) may now be written 
a i Key2(1-n) 4 = @ . hee 1) 


Transformation of equation (4.11) into a Bessel equation can 


be accomplished by the following substitutions. 
a 2, Kk oes Gia 26) 
n 


(4.12b) 


_ 
I 
= 
: 
W 
ae 


Equation (4.11) is then reduced to the form 
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z 
GB. (, - Ys) = 0. (4.13) 
dp p 


where v is representative of 
Se ee (4.13a) 


Verification of the transformation is accomplished in 
Appendix C. 
1. Non-integer Value of v 
For the non-integer values of the parameter v 
between zero end one equation (4.13) has the following solu- 


mon [16] 


— a,J,,(0) = anJ _.(P) (Hee) 


where J, and Des are Bessel functions of the first kind and 
order \v. 1 

Reocaliime 10125), (i.e. 6 = ¥ B(y)), the complete 
solution is determined by 


l-n 


SSS 


ye Metco) 2 crs (0) 1 (4.15) 


The boundary conditions associated with the lateral buckling 


problem are: 


iy = 0) = 0 ee cae 5) (16a) 
= c= a) Pee Ce ee (1) .16b) 


From the first condition and expansion of J_ hie wes COn— 


VY? 
cluded that = O. Therefore (4.15) becomes 
l=-n 
Z e 
o= Y [a,J.(p)} Chel 
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for which differentiation yields 


-l=-n keg 
dé _ l-n 2 2 , dp 
aq?) |"oe eye eae | C18) 


It is additionally known that 


do _ y,1-n 

cca Ky (4.19) 
aaa [17] 

mee = = J,.0p) + 2s / (4.2 

we Re e “v(p) 7 
Mealwuation of - ao thescceond boundary condition with 
y = 1 and p = K/(2-n) yields the solution 

K mk K 
geen ion Oe ‘Bon’ - lean 


The buckling loads may now be determined by solution of 
equation (4.21). This may be accomplished by independently 


calculating each case 
K a 
JI. (——) =0 Cl 2 an) 
and 


RI = 8 CLES ats ) 


US cK) 
l-n vtl ‘2-n 
and determining the initial intersection of the resulting 
Bessel equation curves. For each considered value of n, K 


is found numerically from the argument of J. A computer 


z 

Handbook of Mathematical Functions, National Bureau of 
Standards Applied Mathematics, Series 55, p. 363, June 
1964, 
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solution for this problem is included in the Computer Pro- 
Beem securon of the thesis. The Bessel equation curves 
which show the intersecting roots are included in Appendix D. 
When n=0 the beam design reduces to that of the rectang- 
ular case. The computed solution from equation (4.2la) and 
(4.21b) provides an argument value of 1.05851. Thus, the 


@mitical lead is determined by 


so = 1.05851 (4.22) 
or iM 
oF © 
, 16-936 Kec. 4.23) 
ida We tee 


precisely the value determined previously. The values for 
Seemer solutions are discussed in a later section. 
Oo. Integer Value of v 
When the value of n is equal to one the resulting 
value of v becomes zero and equation (4.13) is reduced to 


the form 
1 wae 
+ —-— —— + = ; tie I 
an = et ee (428) 


Piescouemere Solutiegn ofethis eGquation cam be shown to be 


ot the form [17] 


B = a,J (p) + any (oe) (4,25) 


where 1 is known as a Weber function or Bessel function of 
the second kind. Substituting for ¢, the above equation is 


transtormea to 


Si 





[a,J (po) + any,(e)]. C2 o) 


Mees Solution is valid only if v is of an integer order. 
mms, in the restricted range n is allowedwonly to be equal 
me one, and the y multiplier must be to the zero power. 


Equation (4.24) reduces to 


= 6G) SICES Yanan ye Oey)» Gh 2) 


ie=value of y, (Ky) fworey o= OM ase, x= =) is -© which 


implies that a, is equal to zero. Differentiating $ with 


2 
respect to the remaining term and evaluating at the second 


boundary condition (y = 1, x = 0) requires that 


= UEP (4.28) 


mee the selution to (4.26). This value is found [18] to be 
either the trivial case (i.e. zero) or 3.8317. The latter 
rPeswmiteazmplies a singularity condition with respect to the 
previous non-integer solutions. A graphical representation 
tire. () ot wealwes of mM versus the roots of equation (4.21) 
does, in fact, imply the resulting value converges to zero 
as n approaches one as a limit. Therefore, the trivial case 
is the solution. The beam is apparently unstable for any 


load when in this configuration. 


R. CALCULATION OF AMENDED LOADS 

Extension of the problem is now made in an attempt to 
find the design among all those analyzed that maximizes the 
buciiime iead. Initially ee was @escribed as a fixed height 


at the beam midspan. This center height is now allowed to 
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O v2 4 ‘© 8 1.0 
VALUE OF N 


Pigure {. Convergence of the Root x as N approaches One. 


vary such that each variable design studied will have the 
same volume as the constant rectangular beam. 

This may be accomplished by assuming a constant width 
for all designs, integrating the power law (4.1) and setting 
the resulting value equal to that of the classical rectang- 
ular beam. Using symmetry and integrating the power law 


Trome x=0 to x= * 


L/2 


Den _ 
b ; ho (1 i ax = a Cie edie 


determines the volume of the variable design. For the rec- 
tangular beam the volume for one-half the structure is 


hybd . ; 
“z= Thus, equating the two determines 


ee (n+1) i (4.30) 


It is obvious that the center height must be expanded for 
the variable design in order that the volumes are compatible. 
When De Pootesrrrcbed tO a2 fieed valuewthe critical 


loads may be calculated as described in Section IV-A, i.e. 


a? 





using the intersecting value of (4.2la) and (4.21b). 


2 
Recalling that K = ee the critical load is defined as 
oo 
8 C505 
8(2-n)x,/C7C5 
SS a 


where x represents the intersecting value of the Bessel 


equation curves. The amended load is now found by sub- 


: : O O 
ee lbucing for Ci > C5 and Ao: 


62 cio 
ey io 32) 
US, 
eee =n) 
where 
eo ii 
Cte Ce (4.35) 


The amended load is referred to as P* since its value is 
greater than its unexpanded counterpart. It is left in 
the above form in order that it may be readily compared to 
the solution for the rectangular beam. 

Mm comeucer solugiaon is included in the Computer Program 
Section and provides for determining the intersecting value 
et (4ezta) and (4721p). It additionally calculates the 
critical and amended loads for each value of considered n. 


These results are summarized in Table I below. 


HQ 


TABLE I 


n v vt1 Pew Pp# Payee i 
0 0.2500 1.2500 16 .936=Pu 16 oS ie ao 
al 0.2368 eee 15.618 Ig es up Ol! 
.2 one 22 Iie2e2 14.288 17.145 Le 2 
0.2059 1.2059 12.944 1S oe Th 0.994 
Ges 75 Ig hei i5 les 56.0 IS 2 0.957 
5 0.1667 1 4ai6.67 10.189 15.283 ee 10.2 
6 0.1429 1.1429 8.758 14.013 Ce en 
0.1154 1.1154 (eles 12.347 0.729 
0.0833 aes. 5.656 cE Coaattsel 0.601 
9 0.0455 1.0455 3.794 oe 0.426 
0.98 — 0.0010 1.0010 — be50% i cees 0.059 
ee 0 0 | ~- 0 0 0 


The best design exists between n=0.1 and 0.2. Additional 
calculations indicated that n=0.14 offers the largest buckling 
mee '-=1 7.190) @ror the Variable height formulation. This 
design is shown in Figure 8. Calculation of values for 
P*/Pu indicates that there is little change for the first six 
Values of n. A graph at P* versus n is shown in Figure 9. 
Figures 27-38 in Appendix E illustrate the various designs 


studied in the interval of n=0 to n=l. 


41 





“o6T° LT = rg(WVxe -T)°u = ¥ aseg 
- ret uipeg peguoddns ATduty °9g sansTy 


PLO=N 4S09 -- NOISSG Woes LHOISH 4 lablaun 


GISSHL “NSM “dQ" 1 SN&ld * 1 Loddija 
HONT SLINO OO+300°¢=31N9S~A 
‘HONT SLINM TO+A00*?=31UIS-% 


| 


Ce hod 
co me 


ool 


42° 





P=16.986 2 90% Ry : 
pz mee . weme Oger  _ 
> r) 
10 : 
7 ~o@voam _. = Ne 
= 
© 2 A 6 8 1K 


Figure 9. Amended Load P* Versus Values of N, (case 1). 


fee aNOnON OF THE PROBLEM 


The formulation for design criteria was previously based 


em the assumption 
eel 2x72) Cais 


The case will now be considered for the design based on the 


Conoition 
Reel 2x72) (4.36) 


This condition is physically shown in Figure 10. By inspec- 
tion of the governing equations of the forementioned solution, 


it is evident that the sign change makes no difference in 
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Figure 10. Variable Height Beam Based on h = h,(1+2x/2)". 


mie sOolutTOn ef the problems The idéntieal critical loads 
are again found. However, each design is restricted to the 
same volume of material. This requires that the midspan 
height be decreased while the end heights of the beam are 
mereased. This modification has a significant effect on 
the amended loads. 

Integrating (4.36) and using the same technique described 


in Section IV-B, the center height reduces to 


h Saas ot. eee h (4.37) 
O (gyntl 4 


while the ends expand to 


antl 
me . 2(ntl) ™r (4.38) 
A computer solution for this problem is included at the end 
of the thesis. The results are summarized in Table II. 
The structure is less stable for any comparable condi- 
tion listed in the previous section. Figures 39-48, 


Pepemaix F inelude designs analyzed for this problem. A 
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TABLE II 


n x Per t p# % Decrease 
0 ie 5 1 16.936 16.936 = 
al 102K 7 15.618 ieee 3 tikes Sie 
12 0.99220 14.288 13.215 Pala 
0.95178 12.944 11.508 32.05 
0.90470 11.580 9.891 41.60 
5 0.84907 10.189 8.359 50.64 
6 0.78200 Sue? 58 6.898 Ho, By 
7 0.69840 7.263 5.490 67.58 
8 0.58920 5.656 I Oz (ATES 
9 0.43119 3.794 2.639 gu U1 


plot of the amended load P* versus n is shown in Figure 11. 
Note that all designs reflect a decrease in buckling load 
as compared to Figure 9 and Table I. Also, there is no 
range of designs for which the value P* remains relatively 


ewetionary . 


D. LATERAL BUCKLING OF A CANTILEVER 
A similar analysis [7] for a cantilever beam shows that 


the buckling. loads are determined by 


St (4.39) 


x = ES (4.40) 





a0 





ibs ; 


10 
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Pueure 11. Amended Load P* Versus Values of N (case 2). 


A computer solution was used to calculate the first zero 
Ot -a.meeavive order Bessel function of the first kind, and 
Mee Oeeeioeresult (TO Obtain the critacal and amended loads 
i@mmemecesuruerure. iIhese results are summarized below in 
Table III. The last value is obtained by evaluation of the 
Bessel function of the second kind and use of tables [18]. 
The values of P* versus n are shown in Figure 12. Note in 
particular that all values are greater than the ultimate 
eritical load, Pu, for the rectangular case. 

The best design is determined for the case of n=0.75, 
and is noted as P*=4.918. This shape is shown in Figure 13. 
The range of shapes considered in the formulation are shown 


in Appendix E. 
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epee Ai 


n v Rice an %Z Increase 
0 -.2500 4.013=Pu BeOS —- 
i -.2368 Baca 4.239 5 08 
a -.2222 3.694 4.433 10.47 
3 a9 Bee 4.595 ee 50 
-.1875 3.374 eae eet 
Se) ~.1667 Bee nel ee si2zal eA Oe aS 
Be) -.1429 Sin Olay: 4.885 real are) 
ot ~.1154 geo 4.917 pk ASS 
8 -.0833 eer {Sal Leal 22.48 
9 -.0455 2.568 4.880 PAL 6.0 
me 0 2 405 4.810 19.86 


E. COMPARATIVE ANALYSIS OF SOLUTIONS 

A comparative analysis between the solutions determined 
hor ule Cantilever and simply Supported beams provides some 
interesting insight into the actual physics of the ecient 
The results obtained from the cantilever solution appear to 
be amm@ear Opeimmum shape for the structure. At first sight 
analysis of the simply supported beam does not appear to 
yield the optimum shape, although a better design than 
Piawmaere classical solution 1s obtained. However, a close 
examination of the physics of the problem leads to a predic- 
tion for the possible optimum shape, although it was not 


Mit vacally determined. 
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Figure 12. Amended Load P¥ Versus Values of N 
for Cantilever Ream. 


The value n=0.14 provided a best buckling load of 
pem= 17.199 for the simply supported beam. This corresponds 
to a gain of only 1.55% over that of a constant rectangular 
cross section. For the cantilever, it was mound that 
p¥ = 4.918 (n=0.75) which is a 22.5% gain over fhat on eume 
rectangular section. Physically these results can be 
explained when the lateral bending and twisting moments of 
the two problems are examined. 

Recalling equations (3.4) and (3.6) for the rectangular 


beam yields the following principe axes moments 


5 > (u,-u) ‘ = e =a ae (4.41) 


— 
M" 
| 
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M = = (5 = SG, (4,42) 


Evaluation at the boundary condition, x=0: u,=u and u'=0 
shows that the twisting moment M,=0 and bending moment My 
assumes its maximum value. At the ends of the structure, 
xX = = g?=0 and u=0. This implies that bending moment ee 
and twisting moment Mr is a maximum. 

The combination of these two moments appear to restrict 
the beam from assuming an optimum shape which differs 
maha tmeantiy from the uniform case. In fact, if the 
Beruceures depicted in Appendix E are examined, it may be 
noted that as the beam approaches a pyramidal shape, the 
less stable the structure becomes. Physically, this seems 
EO meply that the vanishing end area offers less and less 
resistance against the twisting moment, and the structure 
will readily buckle. 

it iS apparent and expected that the existing moment 
Genditions govern the optimum shape of the structure. 
Physical intuition would indicate that the maximum amount 
of material be located where the maximum moment existed. 
isis mM and Me were equal and a linear variation could be 
assumed, the resulting optimum shape would actually be 
Genolea recvangular beam. Hor the problem considered, 
MY appears to have a domineering influence on the optimum 
design. If this was not the case, then the solutions 


determined in Section IV-C would have indicated a "best" 


bile lias) load greater than those for the present solution. 
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The actual moment values for M1 and We are not known 
and cannot be determined unless the displacement and twist 
functions u and ¢ are known exactly. Additionally, these 
femctions are coupled to the variable cross sectional area 
as indicated by equations (3.3 a,b,c). Analyzing equations 
eal) and (4.42) and Figure 4 leads to the following inter- 


pretation. The moment representation of M, appears to be 


e 
Darabolic™while that of M, (based on a parabolic ¢) appears 
to be cubic and exhibits an inflection point. These patterns 
are shown in Figure 14. The maximum combined moment must 


then occur between x=0 and x= 5 : 


M = lateral bending moment 


My = twisting moment 


M 
n 





x=0 xe 


Figure 14. Expected Moment Configuration and 
Approximate Optimum Design. 


On the basis of moment evaluation alone, the optimum 
G@esign would apparently be of a form similar to the shape 
shown in Figure 14. Since Mh is expected to be greater 
than Be the design would lump the largest amount of material 


somewhere between x=0 and x= = but closest to x=0. 


By 





Available time has prevented analysis of this expected 
weomcamum' Shape. A Rayleigh Ritz formulation for this 
design was attempted by assuming the displacement and twist 
functions and using a Lagrange interpolating polynomial to 
approximate the design. The results indicated a buckling 
load slightly less than that of the uniform rectangular 
case. 

Because of the low gain in efficiency (1.55%) from the 
best design obtained from the variable height solution, it 
appears that the true optimum shape will not reflect a 
large performance gain over That of the rectangular case. 
In fact, the determined value from the variable height 
solution is expected to be very close to the actual optimum 


design. 


F. VALUES OF N GREATER THAN ONE 

The analyses of the problems discussed have been restric- 
ted to the range for n=l. If n is greater than one the 
actual curvature of the beam design becomes concave. For 
the simply supported beam this physically implies that the 
structure must become less stable than the limiting pyramidal 
shape. This is evident since even less area is distributed 
toward the beam ends. But, it was found earlier that the 
pyramidal shape is actually unstable for any load. By this 
reasoning the problem must be restricted to the range of 
values used to have physical meaning. 

Mathematically, if the value of n is allowed to be 


greater than one, then equation (4.15) changes to the form 
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any [a,J_ (9) + and (Ce) J (4.45) 


ner Which aT muse be Zero. Siew resuering roots introduce 
a singularity condition and the solutions are misrepresen- 
tative of the problem. Consequently, it is deemed necessary 


for the value of n to be restricted to the interval 


feoeee- | to previde meaningful results. 


Sis: 





V. CONCLUSIONS 


For the simply supported beam the assumed power law 
restricts the number of designs which may be analyzed. 
However, a good indication of the problem behavior can be 
determined from the results. The number of designs which 
remain relatively stable do not differ radically from the 
uniform case. As the designs converge to a pyramidal 
shape they rapidly become less stable. 

It appears that the combination of the two moments, 

M and Me govern the design. As the beam approaches a 
pyramidal shape, less area becomes available to resist the 
twisting moment and the beam readily buckles. Thus the 
iMigeracvion of the two moments restrict the beam from 
ameatning a shape which differs significantly from the 
rectangular case. 

The effect of the moment values on beam design is 
Memever Gemonstrated from the cantilever résults. At the 
free end of the beam there is no moment and the design calls 
Poa venye li Guile Material to be located at that position. 
The interaction of the two moments seem to produce a maximum 
slightly away from the fixed end —- otherwise, the design 
for the best buckling load would have indicated the pyra- 
Didatmecmape. Ihe cantilever solution, thus, helps to 


Suwer> the validity of the simply supported beam results. 
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The true optimum shape for the simply supported beam 
has not been determined by the variable height formulation. 
However, the expected optimum shape was approximated from 
woeiaezineg the physies of the probllem. This shape is not 
eeepecued toe offer any significant increasewover the design 


obtained from the variable height formulation. 
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APPENDIX A 


A STRUCTURAL OPTIMIZATION PimebLish 


Po ememple of a structural optimization problem which 
Miwooeacess the usefulness of the variatiemal method of 
Ssectien Il is that of obtaiming the optimum shape of a»euni- 
formly loaded beam. This problem is a "strength" problem 
in contrast to the buckling problem which exhibits no 
eemereun. For the problem considered, bending is in the 


Vertical plane of symmetry (Fig. 15). The height of the 





Figure 15. Variable Width Uniformly Loaded Beam. 


Peoiiets=Gaken as constant, and a uniform load is applied. 
imMewoeavement Of tne Problem 
Wiesdesipn b(x) is to be determined for a variable 
width which provides a maximum strength capacity for the 
structure. Unless some specific constraints are imposed on 
the problem, results may well be meaningless. For instance, 


17 a Minimum volume structure was considered with no other 
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fonGirions than equilibrium and boundary cenmditions to be 
mecioited, then it is entirely reasonable ®o deduce a zero 
volume structure with infinite strength capacity [6]. 

For the strength problem the primary constraint is that 
the material must not fail. Thus, for an elastic structure 
composed of material of finite strength capacity, the 
specific strain energy (v) iain be no greater than a speci- 


fied maximum. This maximum is denoted as w, such una , 


p< (A.1) 


at any point in the structure. The inequality of (A.1) may 
be written as an equality constraint by introducing a scaler 
Peele y(x) to take Wo any existing slack in the function. 


Using the "slack variable" gives 
| <a 
a (A.2) 


For the simply supported beam the moment of inertia with 
respect to the principle bending axis is I = a b(x)h2 or 
T = CA(x). The strain energy per unit length at any posi- 
tion along the beam is known to be 1/2ECAv"* where v repre- 
Pomts tne displacement function. If the strain energy per | 
merc Length is divided by the cross sectional area, A, at 
ir eesmeren, then ’specific strain energy ™per unit aréeam@per 
unit length may be written 
2 


Tee one (A.3) 


An energy constraint equation for the problem is then 


determined, 


BI 





W - (5 ECv"® + y*) = 0. (A. 4) 


Mee cona COnStraint on the preblem restriets the beam to be 
designed from a given volume of material As iis is the 
Meemperimeopric constraint introduced in Section II. Thus 

Q 
(ox) dx =—i (A.5) 
; fe) 


The total potential energy of the system may be written 
a ne 
T = fis ECAV = Paw ax, Cros 


where the first term represents the strain energy due to 
bending and the second term defines the external work of 
the system due to po f10). 

The volume constraint (A.5) and energy constraint (A.4) 
are now accommodated through the formation of the augmented 


minWet ronal T* , 


aaa 12 


T# = f [5ECAv ne 
0 


ae! 2 
— Div A A-Ag CW - SECY 1 Nie Chie 


where aT and do 


Fermerminfma Variation With respect to v, A, and ¥ 


are Lagrangian multipliers. 


yields three Euler equations. 


| 


¥=f1- tt "! = = 
6 T#=0: [ECAv" + ECv",\,] Py 0 (AS 
es 
¥=f)- a '! = 2 
6 ,t OF: ECV Ay 0 (A.9) 
¥=0- = 
6. T=0: 2rny = 0. | (A.10) 


The natural boundary conditions are v(0)=v(L)=0 and 


eo =) )=0. 
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Integrating equation (A.8) twice yields 


O 
1" if) fee 
eecAv' + Aov ) 5 px 7 C,x + C., : Chee aale) 
= 
for which the above boundary conditions determine C ooops 
and C, = 0. Consequently (A.11) becomes 
veer.) = © (x2_¢x) CA ae?) 
O PANG ; ° 
Using equation (A.9) gives 
2A, 7 
va (——) Cle) 
BC 


which physically implies that the deflection curvature is 
Gensrant. 


From equation (A.10) there exist three possibilities 


ili) ho = 0O,y =.0 (A.14) 

ii) +, # 0,y = 0 Ce) 

ee numelay 7 0. - (A.16) 
Assuming y = 0 and examining equation (A.4), it may be 


observed that 
= Ge) (4.17) 
EC ‘ : 


Comparing this with (A.13) implies 


p=) 


cal (Oe ie.) 


Recalling [11] for a simply supported beam under a 
iw toad Fo, the moment at any position X along the beam 


may be represented by 


Mex) = ~2 (x? =. oe) | ON EL) 


5) 





then equation (A.12) may be rewritten 


vi" (Ati, ) = aa M(x). CAL 20) 


But, the moment is also known as M=ECAv"; consequently, 
v"(Ati,) - ee ; Oe ea) 


Prem bine above equation it 18 obvious that A, = 0 and (A.14) 
is satisfied. This condition implies that the energy formu- 
lation actually rejects an energy constraint. This can be 
observed by examining equation (A.8), which should be the 
bauer priwm equation of the system. The eoendition is valid 


oumsy af A, is, in faet, zero. The rejection of an energy 


2 
constraint appears to be inherent to the energy formulation. 
tags formulatwon seeks Towmaximagze thesmaximum stress 

at every HOsmuronex (iter arene energy density) which is 
different than maximizing specific strain energy as was for- 
mulated in the problem. 


Therefore, equation (A.21) becomes 


D 


vex) = see (x ee GO imeele.) 


from which the variable area may be determined 


2 Da (eee tx) 


e O23) 
(8 ~ EC)* 


Lue ) 


Substituting for A(x) into the isoperimetric equation (A.5) 


eG imtesrating, the load Po is found 


po = - — (8UEC)* (A.24) 
QR 
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iiewresult determines the optimum shape of the cross sec- 


Ciemat area, or in terms of the width 


(A.25) 


Since the stress at any section may be represented as 


o = <<, substitution for M and ec (ec ee yields 


5 ee er eee (A.26) 


This represents the maximum allowable stress at any section 
along the beam. It is noted additionally, that the stress 
is independent of any position x along the beam. Thus, the 
stress is maximized at the outer fibers (h./2) for any 
position with respect to the optimum shape. The optimum 
shape is noted to be parabolic and is shown below (Fig.16) 


foarte oumeeses Indicated. = ihesdesien is comparable womche 


uniform strength design listed by Esbach [20]. 





Figure 16. Optimum Shape for Variable Width Beam Design. 
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APPENDEX B 


Breecl OF THE LOAD P ATA DHSTANCE FROM THE CENTROID 


In Section III the assumption was made that the load P 
peould=be located at the centroid of the beam. The effect 
emie Geomeyrical location of the load is significant in 
eme energy formulation analysis. For the lateral buckling 
case, the external work associated with the unstable state 


is represented by 
on J qu’ ($72 = x) dx (220) 


Tres condition results from analyzing a longitudinal element 
ate woeint Din Figure 4. Considering bending of this 
element in the &% plane with the cross section mn assumed 
fimgead, the end of the beam then deseribes an infinitely 


smell arc 
nn Meee” = x) dx (B17) 


mame ool plane, for which the vertical component may be 


represented 
CMe? = xo dx (B.2) 


The external work is P times the integral of this term over 
e@ewarea Considered. 

Now consider the case where the load P is:located at the 
top of the beam (Figure 17). The external work must now be 
determined from lowering of the load by amounts A, and A, 


[19]. | 
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Figure 17. Effect of the Load P Located at Top of Beam. 


A, is simply lowering of the load due to the change in 
Pentroidal height during lateral buckling and is the same 
as (3.20). The component A, represents the lowering due to 


meccineg of the beam throvigh an angle » and can be repre- 


Z 
sented as h/2 (1 - Cos 9) » —— Meelnvic vanveacditiomad Cem 
must be considered if the load is placed away from the beam 


eentroid. This term reduces to ae see 
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Pee en SC 


TRANSFORMATION OF EQUATION (4.11) 








2 
ao,n a, w2y2(l-n) 4 = 9 (4,11) 
dy nf a6 
_ «K Bo (A i2a) 
p oe Day) Y 
l-n 
from (4.12b) 
-l-n l-n 
d@ _ 1l-n 2 2 dp 
where 
d8 _ d8 do 
Spepoocrouuing to get 
=e 2(1-n) 
dd _ i-n 2 2 dp 
dy face 2) ae B 4 mf K do (Cx3)) 
Differentiating again 
5 - 3-n -1-3n 
d“o¢ _ (1-n) (-1-n) 2 2) ade 
a 5 5 4 8 + 2(1l-n)K y ale 
dy 
2 
hi Key 272 (l-n) ad : Cea) 
do 


Peo rmeucing from (4.12b), (C.3) and (C.4) into (4.11) and 


ermraneing in tabular form: 
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d 78 dp : 
do do 
BD 1-—3n ~3-n 
Ke - 2(1=n jiK e né- y e a“s 
Y Y oie — 
— ( ) = do 
z GOCE 2 n Gao 
cal 2 u y dy 
(1-n) _ 
Key? Key? (1 i) 
5 (1-3n) 
9 a(1-n) gee 2 dg 
ie Y m2 a (2-n) Ky do 
0 
-3-n 5 
2 ———— =(l1-n) 
+|(-3-+2- py *% + Key? Be (CaS) 


The operator from equation (C.4) is shown at the right side 
of the table. 


Substituting (4.12a) amd letting 


_ l-n 
ie Sonn CERiG») 


Equation (5) becomes 


2 
ap 1 Gap v 7 
eee Pot] = =) gp =00 Chae) 
dae Oo ap 0° 


Gok 
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APPENDIX D 
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Figure 18. 
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Figure 20. 
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Figure 26. 
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